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SUMMARY 

The observed diffuse galactic y radiation is compared to that predicted 
from galactic cosmic ray interactions with galactic matter and photons, 
assuming that on a broad scale the galactic cosmic rays in the plane are 
correlated with matter oensity. Recent considerations of the qalactic diffuse 
matter distribution, particularly the molecular hydrogen, the galactic photon 
density, and a revised cosmic ray galactic scale height, are included. The 
predictions are compared to the observational y-ray longitude distributions, 
the latitude distribution, and energy spectrum, including the recently 
reported COS-B satellite results, and the new COS-B background estimate. 
Considering the uncertainties, the agreement between the theoretical 
predictions and the y=ray data seems generally reasonable, suggesting that the 
general concepts are likely to be correct. Both the results determined he v '~ 
alone and in conjunction with other work calculating source functions assuming 
only cosmic ray matter contributions Indicate no necessity for a significant 
point source contribution to the diffuse y radiation in the energy range being 
considered (E y > 10 MeV) . The intensity of the highest energy y-rays (Ey > 

300 MeV) could be explained entirely by cosmic ray matter and photon 
interactions, and, in fact, the relatively low intensity of those highest 
energy photons is of some concern In relation to the allowed range of the 
interstellar cosmic ray electron spectrum, If all of the diffuse y radiation 
is to have this origin. Other possible explanations are a larger contribution 
of point sources at low energies, although the few observed y-ray spectra o^ 
point sources would not suggest this, or an enhanced cosmic ray spectrum at 
low energies. The Compton contribution over most regions of the Galaxy is 
calculated to make a 10% to over 20% contribution to the diffuse y radiation 
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depending on latitude and longitude; however, in the inner galaxy the 
contribution is smaller due to the effect of the Compton radiation itself on 
the parent cosmic ray electron spectrum. 


Key Wordiu Gamma Rays, Galactic Structure, Cosmic Rays 

I. INTRODUCTION 

The Y-ray sky is dominated by radiation from the galactic plane, which is 
generally assumed to be the sum of diffuse radiation and point sources. The 
point source contribution would for the most part appear diffuse to the high 
energy Y-ray satellite Instruments that have flown thus far because the 
angular resolution of these instruments for individual photons has been only 
one to a few degrees, or poorer, depending on energy. 

The source of the true diffuse radiation has been assumed to be cosmic 
ray interactions since, assuming cosmic rays pervade the Galaxy, they neces- 
sarily produce high energy Y-rays as they Interact with the interstellar 
matter and photons. The cosmic ray nucleon interactions give rise to y rays 
primarily through the decay of n° mesons, giving a unique spectrum with a 
maximum at approximately 68 MeV. Cosmic ray electrons produce y rays through 
bremsstrahlung, but with a markedly different energy spectral shape, one which 
decreases monotonical ly with energy. Cosmic ray electrons also interact with 
the interstellar starlight, optical and infrared photons, and the blackbody 
radiation through the Compton process. Finally, cosmic ray electrons can 
Interact with magnetic fields giving rise to synchrotron radiation, but this 
process can be shown to be much less Important than the others previously 
mentioned for the galactic diffuse radiation (e.g., Fichte! et al., 1976) and 
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will not be discussed further bare. Should the y radiation due to cosmic ray 
Interactions be dominant for the diffuse galactic emission relative to the 
contribution of point sources, the observations of the diffuse y radiation 
together with other Information should provide a better understanding of the 
general character of our galaxy than would otherwise be possible. 

Substantial work has already been performed on the calculation of the 
source functions for these various y radiations and the Intensity to be 
expected In the vicinity of the solar system. (For a general review see 
Chapter 5 of Flchtel and Trombka, 1981.) There has also been a substantial 
number of attempts to correlate the Y radiation with the matter distribution 
(e.g., Flchtel et a!., 1978; Issa et al., 1981; Stfong et a!., 1982; Arnaud et 
al,, 1982; Lebrun and Paul, 1983; and Riley et al., 1983) with generally 
reasonable results. However, several recent devlopments make a reexamination 
and extension of this work worthwhile. These Include the detailed results of 
high energy galactic Y-radiation obtained with the COS-B satellite (Mayer- 
Hasselwander et al., 1982), further evaluations of the 21 cm radiation in the 
galaxy, and hence the atomic hydrogen density distribution, additional CO line 
observations from which molecular hydrogen column densities are deduced, the 
high photon density estimate for the inner galaxy which affects the Compton 
radiation and the electron spectrum in this region, the current lonqer 
estimate of the galactic cosmic ray lifetime, further evidence support! na the 
galactic arm concept, and Improved theoretical calculations on the nucleon- 
nucleon source function. It should be noted that the difficulty in 
normalizing the molecular hydrogen column density deduced from the CO 
measurements In an absolute manner remains a problem, but one which can at 
least be constrained. 
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An Important assumption of this paper Is that cosmic rays are correlated 
with matter on the scale of arms.,, Because of the complication that this 
assumption Introduces, the reasons for believing It to be correct are reviewed 
again together with the new Information supporting It. It Is worth noting, 
however, that for galactic latitudes where the local contribution may be 
expected to dominate, lb I greater than 10° o' 15°, the cosmic ray density as a 
function of the local galactic position may not vary much, For this case, 
since the scale height of the cosmic rays is expected to be large compared to 
that of matter, a good approximation for the cosmic ray, matter interaction 

t 

contribution to the Y-ray diffuse radiation Is probably obtained by using a 
constant cosmic ray density, which allows the direct use of atomic and 
molecular hydrogen column densities. If the point source contribution is 
small arid If account is taken of the Compton contribution, It should be 
possible to obtain a good agreement using the matter column densities directly 
as shown by Strong et al. (1982). It should also be possible to use this 
simplified approach successfully at Intermediate longitudes (~ 60° to ~ 100° 
and ~ 250° to ~ 280°), where regions which are at galactic radii similar to 
the Earth are predominantly being viewed, as shown, for example, by Arnaud et 
al. (1982) and Lebrun et al. (1983) for the (60° < 1 < 100°) region. The 
developments related to the galactic matter distribution and Y-ray production 
will be considered In the next section and Incorporated Into the general Y-ray 
production calculation. The predictions of this work are compared to the 
recently published high energy Y-ray results in Section III, with the 
conclusions summarized in Section IV. 
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II. DIFFUSE GALACTIC GAMMA RAY PRODUCTION 

i (a) Galactic Matter Distribution 
With regard to the matter, the relevant concern Is the qalactlc diffuse 
matter In the form of atoms, molecules, 1ont£ and dust. The latter two are 
believed to be minor constituents and, hence, unimportant for y-ray production 
through cosmic ray Interactions. Hydrogen Is the primary component of both 
the atomic and molecular matter. Helium and heavy nuclei add about 05% more 
to the v-ray production. It Is assumed these latter nuclei have a distri- 
bution In the galaxy similar to hydrogen, although little Is known about them. 
Both atomic and molecular hydrogen are known to be confined to a narrow disk 
with the molecular hydrogen distribution generally having a smaller scale 
height (e.g., Gorddn and Burton, 1976; Solomon and Sanders, 1980)* 

If It were true that the cosmic ray density were constant throughout the 
galaxy, It would only be necessary to know the column density of the hydrogen 
In order to calculate the diffuse galactic y-ray emission (See, for example, 
Flchtel and Knlffen, 1974). However, If the cosmic ray density Is variable, 
the product of the cosmic ray density and the matter density must be inte- 
grated over the line of sight in the galaxy, and hence, the matter distri- 
bution In the galaxy must be deduced. This point will be discussed further in 
Section 11(c). It Is worth noting now, however, that for the local region of 
the galactic disk, represented by |b| *,10°, where the cosmic ray density In the 
plane varies slowly, and for regions where the Compton component variations 
are not too Important, the galactic y-ray s themselves represent an indication 
of the galactic column density as a function of direction. Correlation 
studies of this type have been performed by Lebrun et al. (1982), Lebrun and 
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Paul (1982), Strong and Wolfendale (1981), and Strong et a]. (1982), and their 
relationship to this work Is discussed later. 

Consider first the neutral atomic hydrogen. Its density as revealed by 
the 21 cm emission remains somewhat uncertain In the Inner galactic regions 
because of uncertainty In the absorption correction. Recent work (e.g., 

Dickey et al., 1982; Thaddeus, 1982) suggests that the absorption had 
previously been somewhat underestimated and that the density In the region of 
3 to 5 kps from the galactic center Is probably greater than previously 
estimated perhaps by a factor of 1 1/2. In this work, the atomic hydrogen 
density distribution of Gordon and Burton (1976) as a function of radius from 
the galactic center was usea, but modified so that the atomic hydrogen density 
in the Innermost region was increased by a factor of 1.5, and the closer 
densities were increased less In accordance with the amount of intervening 
matter. In earlier work (Knl f fen, Fichte! , and Thompson, 1977), a scale 
height of 0.12 kpcs had been used In the Inner galaxy and a value gradually 
increasing from 0.12 kpcs in the outer galaxy. Recent work by Lockman (1982) 
has shown, however, that the effective scale height is about 3/2 times larger 
then previously believed because a relatively faint component of HI has been 
overlooked. Hence, the scale height used in this work is 0.18 kpcs for the 
galactic radius, R Ga -| , less than that of the Sun (10.0 kpcs) and [0.18 + 0.023 
x ( R Gal “ 10.0)3 kpcs for Rggi > 10.0, with the increase beyond the solar 
circle being based on the work of Baker and Burton (1975). It Is now believed 
that the scale height In the outer galaxy Increases more rapidly than this 
(Kulkarni, Blitz, and Heiles, 1982), but the surface density, which is the 
relevant parameter whervthe total galactic plane contribution is considered, 
is still believed to be similar to earlier estimates. However, this larger 


scale height would Imply a somewhat broader latitude distribution In the outer 
galaxy. The latter authors have also shown that the galactic disk extends to 
30 kpc although, by that distance, the surface density has become quite small, 
of the order of 0.1 M 0 pc” 2 compared to 6 M 0 pc” 2 near the solar system. The 
density distribution used in this work Is also modulated for the qalactic arms 
1r« a tanner to be described below. 

The density distribution of molecular hydrogen Is measured less directly. 
At present, the best estimate is obtained through the observations of the 
2,6 mm spectral line of * 2 C0, from which the distribution of cold Interstellar 
matter is inferred. The nature of the Interpretation of these measurements 
makes the derived molecular hydrogen density distribution less certain than 
that of the atomic hydrogen. The average galactic radial distributions of 
molecular and atomic hydrogen show clearly that the molecular hydrogen to 
atomic hydrogen ratio is larger in the Inner galaxy than It Is in the outer 
galaxy even If the absolute intensity of molecular hydrogen Is still quite 
uncertain. The basic distribution of the density as a function of distance 
from the galactic center was taken from Robinson et al. (3.983) and Gordon and 
Burton (1976) with the molecular hydrogen density normalization treated as an 
adjustable parameter. (It is interesting to note that the CO observations 
indicated that the great majority of the molecular hydrogen Is In clouds. The 
work of Solomon and Sanders (1980) has, In fact, suggested that the Inter- 
stellar medium is dominated by massive cloud complexes.) 

Although the translation of the observations into a galactic spatial 
distribution is difficult, on a broad scale the density profile is reasonably 
well accepted. Even though there is no general agreement on details of arm 
structure, a general spiral pattern does appear to emerge. In addition to the 
21 cm data, the distributions of continuum radiation (Landecker and 



Wleleblnskl, 1970; Price, 1974), Y radiation (Blgnami, et at., 1975), Hill 
regions (Georgelln and George! In, 1976), supernova remnants (Clark and 
Caswell, 1976), pulsars (Selradakls, 1976), and Infrared emission (Hayakawa et 
al., 1976) are all consistent with the existence of spiral structure In the 
galaxy* Until recently, It had not been clear whether molecular clouds were 
associated with spiral structure. However, now on the basis of a high sample 
survey and observations In both the first and second quadrants of the galactic 
plane, Cohen et al. (1980) have reported the existence of the molecular 
counterparts of the five classical 21 cm spiral rams segments In these 
quadrants, namely the Perseus arm, the Local arm, the Sagftarlus arm, the 
Scutum arm, and the 4 kpc arm. Kutner and Mead (1981) have even Identified 
arms through CO measurements In the outer galaxy. The specific spiral pattern 
that will be used here is that of Georgelln and Georgelln (1976). In regard 
to the particular choice, the spiral structure model recently developed by 
Robinson et al. (1983) based on a current well -sampled CO survey by these same 
authors shows “excellent agreement" with the Georgelln and Georgelln model. A 
five hundred parsec width Is adopted for the arms. The excess of material In 
the arms Is taken to be twice the local average density of matter not In the 
arms, unless the distance between the arms Is less than the arm’s width in 
which case It is proportionally smaller, based on recent considerations 
(Lockman, 1982, and Kulkarni, Blitz, and Heiles, 1982). In either case, the 
total matter Is made to be consistent with the estimated column density, 
although this is practical only on a broad scale and not on a fine scale (of 
clouds, for example). 

It Is realized that the unfolding of the radio observations to obtain a 
spiral arm, or arm segment, matter distribution pattern fs necessarily 
somewhat uncertain. However* If the theory is to incorporate properly the 
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correlation of the cosmic ray density with the matter density on the scale of 
arms, It Is necessary to have the density dl stlbutl on and not just the column 
density. A reasonably accurate calculation of the T-ray Intensity does not 
require an exact galactic picture of the matter distribution because a good 
approximation will result If the mass Is In approximately the right place, and 
reasonable care has been taken to verify that the arm width and mass ratios 
do, In fact, lead to column densities that are In agreement with the column 
densities deduced from observations on the average. Some fine details on the 
scale of clouds are lost. 

(b) Galactic Photon Distribution 

For the photon distributions, Knlffen and Fichte! (1981), using 
results of Boissd et al. (1982) on the Infrarsd volume emlsslvlty and a model 
of Bahcall and Sonelra (1980) for the starlight distribution, obtained photon 
densities and, hence, a source function for the Compton emission as a function 
of position In the galaxy* These will be used here. 

(c) Galactic Cosmic Ray Distribution 

With regard to the cosmic ray distribution In the galaxy (see 
particularly Knlffen and Flchtel, 1981, and Flchtel et al., 1976), It will be 
assumed that the nucleonic cosmic ray composition and energy spectrum remain 
unchanged throughout, the galaxy and that the electron spectrum changes only In 
a second’ order manner as the density changes, except at very high energies 
principally In the inner galaxy. The latter point will be discussed later in 
Section 11(d), The cosmic ray density In the plane will be assumed to be, 
proportional to the matter density on the scale of arms and, perpendicular to 
the plane, to have a Gaussian distribution with a scale height of 0.6 kpc. 

This value Is based on the radio continuum measurements of Cane (1977) and the 
assumption that the galactic magnetic fields energy density and the cosmic ray 
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energy density have the same scale height. This scale height for the cosmic 
rays Is somewhat less than that use$ previously, and the primary effect 1$ 
some reduction In the Compton contribution. 

Since the assumption that the cosmic ray density Is correlated with 
the matter density substantially complicates the calculation, It Is worth 
reviewing why it Is believed to be the case. The galactic cosmic ray pressure 
locally Is about equal to the magnetic field and kinetic motion of matter 
pressures, and together they are as large as can be held by the local galactic 
matter. These conditions suggest that the cosmic ray density is as great 
locally as the galactic matter will allow. Further, the cosmic ray age deter- 
mination suggests that this situation Is the result of plentiful sources and 
leakage, not just chance accumulation to the maximum over time. Hence, 
excluding the possibility that the local conditions are anomalous, the most* '' 
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natural assumption is that the cosmic ray pressure is as great as it can be 
throughout the galaxy except possibly in the outer galaxy where sources or 
regions of further acceleration may be rare, (For a further discussion, see 
Chapter 5 of Fichte! and Trombka, 1981.) The assumption that the density 
varies on the scale of the arms is based not only on the natural scale of the 
arms, but on the scale height of cosmic ray electrons perpendicular to the 
plane ~ 600 kpcs and the theoretically suggested mean diffusion length in the 
plane (a few to several tenths of a kiloparsecs). Support for this assumption 
is obtained from the recent work showing that the cosmic ray electron 
Intensity within the spiral arms is about a factor of 2 higher than between 
the arms (Webber, 1983) \ 

It is ultimately the total gravitational mass that is relevant in 
considering the galactic attractive force needed to balance the expansive 
pressures of the cosmic ray gas, the magnetic fields, and the kinetic motion 
of matter, e.g., Parker (1966), and there is much more mass in the stars than 
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In the diffuse matter* The assumption that the cosmic ray surface density is 
proportional to the diffuse mass surface density In particular was made 
because the distribution of the latter In our galaxy Is better known. The 
observations of other galaxies Indicate that the two populations are similar 
at least for Rg^ < 15 kpc, and the galactic magnetic fields which control the 
cosmic ray motion are generally believed to be correlated with the diffuse 
mass. 

The combination of the simp'. . -ylng assumption of step function matter 
arms and a similar cosmic ray distribution Is recognized to be unrealistically 
sharp and leads to enhancements In the predicted Y-radlatlon from the direction 
of arms which are too sharply defined. However, the refinement of a smooth 
rise and fall would be difficult to implement, and this seems unwarranted at 
this time both on the basis of the limitations In the Y-ray data and the lack 
of knowledge of tfti* setter arm profile. 

(d) tes# Ray Source Function and Calculation of Predicted Intensities 
The detailed calculations associated with the production of 
energetic y rays through cosmic ray nucleons interacting with Interstellar matter 
including all the primary cosmic ray and Interstellar matter components, all 
the secondaries and their decay products, the angular distribution, and the 
energy spectrum are very detailed and lengthly. These calculations have, 
however, been performed. Following the original work of Cavallo and Gould 
(1971) and Stecker (1971), Badhwar and Stephens (1977), Stephens and Badhwar 
(1981), and Morris (1982 and 1983) have used the substantial recent hloh 
energy physics experimental work to estimate the Y-ray production energy 
spectrum for cosmic rays Interacting with Interstellar matter. The spectral 
shapes calculated by Stecker (1973), Badhwar and Stephens (1977), and Morris 
(1982) with the corrections at high energies (Morris, 1983) are In fact very 
similar. The Morris (1983) work was used here, and the relevant values are 
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given In Table I at the end of this section. 

The cosmic ray electron, matter Y-ray production can be calculated using 
the bremsstrahlung cross-section formulas of Koch and Motz (1959). The 
predicted radiation In the region below about IQ 2 MeV is uncertain even 
locally In our galaxy because the Interstellar cosmic ray electron spectrum Is 
not well known at low energies where the electron spectrum observed near the 
Earth has undergone significant solar modulation. 

To consider the electron spectrum, It is necessary to look ahead to some 
of the y-ray results. At present, the Y-ray spectral observations, particu- 
larly of COS-B, represent a concern Independent of questions of the details of 
the matter density and distribution and of the variation of the cosmic ray 
density. Prior to the recently reported results of COS-B in the energy region 
above 300 MeV, the galactic Y-ray spectrum had been consistent with a galactic 
cosmic ray population consisting of the sum of cosmic ray nucleon spectrum and 
the electron energy spectrum observed at high energies ( e. g., Protheroe, 1982), 
approximately 3.9 x 10” 3 E~ 2 * 3 electrons cm” 2 sr“* s~* MeV“* changing to a 
spectrum of the same form, but with a power law Index of 2,1 to 2.3 at (Z to 3) 
x 10 3 MeV (e.g., Knlffen and Fichte! , 1981, and Webber, 1982) In agreement 
with the Interstellar electron spectral range set by radio observations 
(Cummings et al., 1973) and with current concepts of solar modulation (See 
also Webber, Simpson, and Cane, 1980). This combined spectrum changes little 
In shape with position in the galaxy whether one assumes a cosmic ray density 
proportional to matter or a constant cosmic ray density If the life-time is 
similar to what it is locally. (For a discussion of this point, see Fichte! 
et al., 1976.) If this spectrum is used, although there is good agreement in 
shape in the medium energy range (a few to 35 MeV) and over the 35 MeV to 
300 MeV range, where results of SAS-2 and COS-B are available (e.g., Bertsch 
and Knlffen, 1983), the predicted intensity above 300 MeV exceeds that 
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observed by COS-B, Lebrun and Paul (1983) have deduced a spectrum of the 
form E~ 2 * 7 for the energy range below 1 GeV from the Y-ray data alone. A 
cosmic ray differential electron spectrum with the form E’ 2,ft and a normali- 
zation below (somewhat more than a factor of two) that of the observed cosmic 
rays could be used with the cosmic ray nucleon spectrum to obtain agreement 
with the medium and high energy spectrum (Lavlgne, 1982); however, in addition 
to Its being inconsistent with the observed ratio of high energy cosmic ray 
nucleons to electrons, It Is only marginally consistent with the electron 

spectrum deduced from radio data especially when one considers shape. 

# 

For this work, the observed high energy cosmic ray electron spectrum of 3.9 
x 10 3 E” 2 *® electrons cm" 2 sr” 1 s" 1 MeV -i has been retained above 4000 MeV, 
and a spectrum of 324 E" 2 *® electrons cm” 2 sr"* $”* MeV“*, confined to agree 
with the spectral range deduced from radio measurements, has been used below 
that energy. As will be seen, the resulting spectral agreement with obser- 
vation is fair, but not entirely satisfactory. It does not appear possible to 
obtain better agreement If one retains the constraints of not reduclnq the 
high energy electron spectrum below the locally observed value and staying 
within the bounds set by radio observations. (For a further discussion of the 
electron spectrum, see Webber, 1983) This general problem will be addressed 
further In the discussion section, particularly with respect to possible 
alternative explanations. However, It should be noted that If the threshold 
of the high energy range of 300 to 5000 MeV of COS-B were In fact only (30 to 
50) MeV higher than quoted, the difficulty would not exist since the 
theoretically predicted intensity would then be (11 to 17) % lower. There is 
no reason known to the authors, however, to suspect such a correction is 
appropriate. 

It should be noted that to obtain agreement at low Y-ray energies (Sacher 
and Schoenfelder, 1982) there has to be an Increase In slope of the electron 
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spectrum in the region below about 30 to 50 MeV, but this increase is permitted 
by the radio measurements, and it is not unreasonable that such a low energy 
electron component should exist, for example, as £ result of stellar flares. 

The calculations associated with the production of Compton y rays have 
been performed in some detail for cases of astrophysical Interest by Ginzburg 
and Syrovatskii (1965). Cosmic ray electrons Interact with galactic starlight 
photons, for which the optical and Infrared ranges are the Important ones, 
and with the universal blackbody radiation. The source functions of these 
interactions are much smaller in the galactic plane than that for-bremsstrahlung. 
However, the total contribution to the galactic y radiation is significant 
because the cosmic ray and stellar photon scale height above the galactic 
plane are greater than those of the matter. 

In early work, the bremsstrahlung /and Compton spectra were calculated 
from the electron spectrum in a manner described by Fichtel et al. (1976), and 
it was shown that within the expected range of average mass densities over the 
galaxy, although the cosmic ray intensity might vary within a limited range 
according to the model, the spectral shape did not change in a significant 
manner. However, there are two new considerations, one related to galactic 
photon densities and the other to cosmic ray lifetime, which now enter the 
picture, and, although they do not effect the bremsstrahlung predictions, they 
do effect the predicted Compton radiation in the inner galaxy. Knlffen and 
Fichtel (1981) calculated the expected photon densities through the galaxies 
based on existing observations and found that they were much higher than 
anticipated. Particularly in the optical and infrared regions of the spectrum 
they ranged from 3 or 4 to 10 or more times their local value in the inner 
Galaxy (R < 5kpc). The significance of these estimated photon densities lies 
In the effect they have on the very high energy electron spectrum and hence 
the Compton radiation. In the y-ray energy range of interest herr there is ’ 
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no significant affect for bremsstrahlung which for cosmic ray like spectra 
arises very promi nantely from electrons with energies In the range of one 
decade Immediately above the Y-ray energy. The Compton radiation typically 
comes from much higher energy electrons. 

To approach this point more quantitively, consider the quantity 

t 

f e e \ 


for synchrotron radiation, Compton radiation, bremsstrahlung, and ionization, 

where is the electron lifetime, and E e Is the electron energy. If this 

quantity approaches minus one for any region In the galaxy In the energy range 

of Interest (for the purposes of this discussion E fi > 50 MeV), the assumption 

of a nearly constant spectrum is not valid. With regard to bremsstrahlung and 

ionization there is no concern in this higher energy range (although for lower 

energies there Is). Using the work of Flchtel et al. (1976), but a t value 

© 

of 1.5 x 10^ years rather than the lower value they used to reflect the fact 

that current thinking places the lifetime of cosmic rays in the (1 to 2) x 10? 

years range gives: 


T e dE e -5 B i 2 

synch BTnocaTT E e- 


( 1 ) 


where B x cj oca T) is assumed to be 3 x 10~ 6 gauss, and E e is in MeV, and 


T e dE e 

r r r — 

L e Comp. 


>4.4 x 10’ 5 


M ph 

y ph(local ) 



( 2 ) 


where B x is the perpendicul ar component of the magnetic f1eld,E e Is in MeV, 
is the photon density, and (] C/Ca u is estimated to be 1.16 cm -3 for 
infrared, optical, and blackbody combined. Other contributions to the Dhot.on 
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density are assumed to be small. Assuming the local electron energy spectrum 
to be in equilibrium, equation (2) shows clearly that in the inner galaxy 
where 

Vi >> M ph( local)* ^ or E e ~ 2 x 10 ^ MeV cosm * c ray electron intensity 
will be dramatically reduced, and for the energy range from approximately 
0.4 x IQ 4 MeV to 2 x 10 4 MeV the effect will be strong. Assuming B^ 2 is 
generally proportional to the mass, Compton radiation losses will dominate in 
the inner galaxy and synchrotron losses may be neglected for the.jP1r$t order 
consideration relevant here. The Compton radiation will Itself be seen to be 
small compared to that from matter Interactions. 

For the energy ranges of particular Interest here. (70 MeV < E y < 160 MeV), 
(150 MeV .< £ f < 300 MeV), (E y > 100 MeV), and (300 MeV < E y < 5000 MeV), the 
typical parent electron energies (although, in fact, there is a broad range in 
each size) are in the range from 0.4 x 10 4 MeV to 1.1 x 10 4 MeV for starlight, 
1.6 x 10 4 to4.4 x 10 4 MeV for the infrared region, and 2.0 x 10 5 to 5.6 x lO* 5 
for the blackbody radiation. As a result, the electron spectrum will be 
sufficiently depressed in the inner galaxy in the relevant energy ranae such 
that there will be essentially no black black body and little infrared Compton 
radiation in this region. Even the starlight Compton radiation will be 
reduced to some degree. In the ir»ner galaxy, were it not for this effect, the 
Infrared Compton radiation would be twice or more that of the optical, and the 
blackbody Compton radiation would be about two-thirds that of the optical 
region. The net result is that the Compton contribution is quite small. This 
effect causes the center to anticenter diffuse galactic Y-ray emission ratio 
to be less than it would be If it did not exist. 

The source functions obtained in this work are given in Table I. It is 
interesting to compare these values with those of other authors making direct 
comparisons of the Y~ray intensities to the matter column densities. Strong 
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at al. (1982) obtain (1.40, 0.53, and 0.59) x 10" 26 y" 1 (H Nuc)" 1 s** 1 ster** 1 
for the three energy Intervals (70-150) MeV, (150-300) MeV, and (300-5000) MeV 
for a sum of 2.52 x 10“ 26 Y (H Nuc)" 1 s" 1 ster" 1 , for 10° < |b| < 20°, 

Issa et al. (1981 obtain 2.2 x 10" 26 y (H Nuc)" 1 s" 1 ster" 1 for E y > 100 MeV. 
For the sake of comparison, the numbers In Table 1 must be divided by 4 * to 
obtain q/4 tr. Doing so, one obtains (0.88, 0.55, 0.53, 1.96, and 1.52) 
x 10" 26 y (H Nuc)" 1 s" 1 ster" 1 for (70 <E< 150), (150 <E< 300), (300 <E< 5000) 
(70 <E< 5000) and (E > 100) MeV, respectively. These latter values are 
generally lower as they should be, since they refer only to cosmic ray matter 
Interactions to which Compton Y radiations and point source contributions are 
added, whereas the former values have knowingly Ignored these contributions 
and attributed all the radiation to matter. For the (70 to 500) MeV region, 
the ratios of the above number are 0.78 and 0.92, and for the (Ey > 100) MeV 
comparison, the ratio of the above numbers is 0.69. The average 0.80 would 
suggest a Compton and point contribution of typically 20%. For the highest 
energy Interval, the two approaches lead to almost the same source function 
suggesting a minor contribution from point sources and Compton radiation. 

In this work, the Compton contribution averaged over |b| < 10° In these 
energy ranges varies from about 6% to about 14% of the radiation due to cosmic 
ray, matter interactions depending on the energy range and longitude, 
suggesting that essentially the entire galactic y-ray diffuse radiation In 
this broad energy range can be explained by the sum of cosmic ray nucleon- 
nucleon Interactions, bremsstrahlung, and Compton radiation, with little 
requirement for the addition of point sources. The ratios of the Compton 
radiation to that from the sum of the cosmic ray, matter Interactions for 
several typical directions and energy intervals are given in Table II. 
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TABLE 1 


Energy Range (MeV) 

70-150 

150-300 

300-5000 

> 100 

Nucleon-Nucleon 

3.4 n B r 

4.5 n B r 

5.6 n^r 

\ 

12.3 nn r 

Bremsstrahlung 

7.9 n H r 

2.4 

1.2 n M r 

6.8 n B r 

Cosmic Ray, Matter 

11.1 n B r 

6.9 n H r 

6.7 fi H r 

19.1 n H r 

Compton Optical 

0.26 Myj S r 

0.12 Mvi s r 

0.13 My<| §r 

0.38 y vis r 

Compton Infrared 

*• 0.44 u IR r 

0.20 Ujpr 

0.22 Mjgr ^ 

0.65 Pjpr 

Compton Blackbody 

0.14 u BB r 

0.07 u BB r 

0.07 u BB r 

0,21 u BB r 


Gamma ray source functions in units of 10“ 26 Y-rays cm“ 3 s" 1 for the energy range 
indicated. "n H " represents the number of hydrogen nuclei per cm 3 either in atomic 
or molecular form. The matter source functions all Include a correction for helium 
and heavier nuclei as described in the text. "Uyj S ", "ujg", and "u BB " are the 
photon densities for the visible, infrared, and blackbody ranaes, respectively, "r" 
is the ratio of the cosmic ray spectra to their local value If the spectra are 
unchanged. For the inner galaxy, the situation for Compton radiation is comnlex as 
explained in the text because of the high rate of energy loss of the parent 
electrons. The values of the 'Vs 1 ' as a function of position in the galaxy are 
given in Kniffen and Fichtel (1981). 

Given the source functions, the Intensities in any direction are then 
calculated in a manner described, for example, by Fichtel and Trombka (1981) with 
consideration of the angular resolution of the instrument being taking into account 
where appropriate. 
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III. GAMMA RAY RESULTS AND THEIR INTERPRETATION 

The predicted T-ray Intensities are compared to the SAS-2 and COS-B longitude 
distributions in Figures 1 and 2. It was noted earlier In this article that 
molecular hydrogen density normalization was left as an adjustable parameter. 

In earlier work (Knlffen, Fichte! , and Thompson and Kniffen and Fichte! , 1981) a 
normalization of 0.6 was used relative to Gordon and Burton (1976). This 
normalization Is consistent with the Independent analysis of radia data by Blitz and 
Shu (1980), and corresponds to a value of 1.3 x 10 20 mol. cm” 2 K” 1 
km"* s *. This value seemed again here to be about the best considering the various 
constraints set by the data in spite of the several new features. There were two 
competing affects in the new consideration. A higher value would have been 
Indicated by the smaller arm to inner arm matter contrast (which has the effect of 
reducing the Intensity because of the assumed correlation of the cosmic ray surface 
density with the matter surface density leading to a Y-ray Intensity proportional to 
the square of the matter density) and the smaller Compton contribution due to the 
smaller scale height for cosmic rays assumed here relative to the earlier work. 
However, the relatively small diffuse galactic center Y-ray Intensity reported by 
Mayer-Hasselwandeer et al. (1982) above 300 MeV dictates that the molecular hydrogen 
normalization factor be kept smaller. 
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TABLE II 


ENERGY INTERVAL 

1 

b 

(70-150) 

MeV 

O 

o 

o 

O 

II 

ii 

o 

O 

o 

O 

fH 

II 

ii 

90° 

o 

o 

II 

K 

90'* 

o 

O 

H 

II 

ii 

o 

O 
f'- 
— i 

O 

O 

II 

ii 

O 

O 

N 

H 

o 

O 

•H 

II 

ii 

290° 

O 

O 

II 

ii 

290° 

o 

O 

iH 

(150-300) 

MeV 

O 

O 

o 

O 

it 

ii 

o 

O 

o 

O 

ii 

« 

o 

O 

CM 

o 

O 

M 

ii 

o 

O 

CM 

i— » 

o 

o 

ii 

ii 

1— ■ 

O 

o 

O 

o 

m 

ii 

o 

O 

r- 

10° 

ii 

ii 

290° 

O 

o 

ii 

ii 

290° 

►-* 

o 

o 

(300-5000) 

MeV 

O 

O 

o 

o 

ii 

ii 

o 

O 

o 

O 

tH 

ii 

ii 

LO 

o 

o 

o 

O 

ii 

ii 

90° 

o 

O 

i—i 

ii 

it 

O 

O 

I s - 

O 

o 

n 

ii 

M 

o 

o 

10° 

ii 

ii 

290° 

o 

o 

ii 

ii 

290° 

o 

O 

i—4 


Ratio of Compton Radiation to that 
Produced by Cosmic Ray, Matter 
Interactions 

.04 

.13 

.06 

.14 _ 

.07 

.10 

.05 

.23 

.03 

.10 

.05 

.12 

.06 

.08 

.04 

.20 

.03 

.09 

.05 

.12 

.06 

.08 

.04 

.20 
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Regarding specifically Figure 2 and the energy spectrum shown in Figure 3, 
there are two general comments. First, the Caravane collaboration has 
recently reanalyzed the combined Instrumental and as trophy si cal isotropic Y- 
ray background (Strong, 1982). Whereas the background intensities for the 
(150-300) MeV and the (300-5000) MeV energy Intervals are essentially 
unchanged from the values given by Mayer-Hasselwander et al. (1982), the (70- 
150) MeV background intensity Is now estimated to be 4 x 10" 5 photons cnr 2 s” 2 
ster“* rather than 3.2 x 10"® photons cm~ 2 s“*ster~*. This chaoge«.has been 
Incorporated in Figures 2 and 4 by Introducing a new "zero base". Second, In 
general, considering the difficulties and pioneering nature of the 
experiments, the agreement between the SAS-2 and COS-8 data is remarkably qood 
in terms of (general intensity level, energy spectra, and relative 
distribution, as seen for example in Figs. 1, 2, and 3. 

Considering the uncertainty in the point source contribution and the mass 
distribution, the agreement between the data and the predicted curves seems 
reasonably good especially when the sources noted In Figure 1 are taken into 
account, except for concern about the energy spectrum which was alluded to 
earlier. 

There are some specific features to be noted in Figures 1 and 2. Notice 
that the edges of the Sagittarius and Crux arms at about 55° and 310° 
respectively mark the beginning of the higher Intensity associated with the 
central region of that galaxy, and that further steps near 35° and 330° mark 
the edges of the Scutum and Norma arms. There appear to be Increases at J5° 
and 285° associated with the local arm and the Carina arm respectively. The 
expected Increase at 265° for the local arm is masked by the large Increase 
due to the Vela' pulsar. / 
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ORIGINAL PAGE 19 

OF POOR QUALITY 
tne model 


have been calculated 


The latltu' distributions resulting from tf 
taking into accot \ the COS-B Instrument response* Figure 4 shows the 
predicted latitude distributions at longitudes near the galactic center, the 
anticenter and at two Intermediate longitudes plotted together with the COS-B 
observations (Mayer-Hasselwander et al., X9$2) . The dashed line In the 70-150 
MeV energy range again Indicates the adjustd zero level on top of which the 
model predictions are plotted to account for the new higher background 
estimate (Strong, et al., 1982). In comparing with the observations, It 
should be remembered that the model makes no attempt to Include contributions 
due to local clouds. The fits to the galactic center observations (1 * 350° 
to 10°) and In the Intermediate range In the first quadrant 40° > 1 > 70° are 
very good at medium and high energies. At (1=275° to 320°), the latitude 
dependence has the right shape, again Indicating the radial distribution of 
the matter In the model Is consistent with the Y-ray observations; there Is a 
slight displacement from zero latitude at all energies In the observations 
(the H hat brim" effect) which the model does not attempt to reproduce. The 
anticenter result Is also reasonable In these energy intervals considering the 
uncertainty of the Y-ray measurements and the gas distribution In that 
direction. In general, these results give confidence that the spatial 
distribution of the emission calculated from the model must be approximately 
correct. The predicted Intensity In the 70-150 MeV range Is somewhat low by 
an amount which appears to be Independent of latitude and longitude. Since 
the shape of the distribution Is In generally good agreement with the data In 
this energy range, except for this constant displacement, the presumption Is 


that the difference is less likely to be due to an underestimated cosmic ray 
electron intensity In the 70-150 MeV range than to some unexplained background 
or nearly Isotrlpic component. Otherwise, In consideration of the 
uncertainties Involved, the overall agreement seems good. 



It should also be mentioned that In the concept being presented here the 
arms on the far side of the galaxy make an Important contribution for small 
(|b| ;C 0,4°) galactic latitudes. This feature results from the high energy 
7 -rays being essentially unattenuated as they pass through the galaxy to the 
Earth and the far side arm matter being largely concentrated In this small 
latitude Interval. With future high resolution Y-ray measurements, these back 
side arms should appear as a narrow ridge superimposed on the broader ridge of 
the near side arms. If the mass density and photon density In the far side 

** 'XSfWm 

arms are similar to the near side arms, although the total Intensity over the 
arm width will be lower by the ratio of the distances (not the distance 
squared because they are arms, not points), the peak Intensities will be 
approximately the same since the distance factor Is canceled by the area 
factor for a uniform density region In a given solid angle. Regions In 
longitude, not directly towards the center, but before the first arm tangents, 
e.g., (5° x, 1 <• 15°) would be ones In which to look for this effect. It 
should also be possible to Identify tentatively very large far side molecular 
clouds If the majority of the molecular hydrogen Is In large clouds. 

It was noted earlier In section lid that the observed Y-ray energy 
spectrum presented a difficulty In terms of the cosmic ray spectra used here 
also having to agree with both the observed cosmic ray Information and 
constraints placed on the cosmic ray electron spectrum by the radio data. It 
Is seen that these constraints do permit fair agreement with the Y-ray data, 
but do not permit quite as steep a Y-ray spectrum as reported by Mayer- 
Hasselwander et al (1982) based on the COS-B data. 

In section IIA, It was mentioned that for certain directions the cosmic 
ray density might be expected to be reasonably constant over that portion 
of the llne-of-slght Integration making a significant contribution to the 
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r radiation and for these regions it would be possible to simplify greatly the 
calculation by assuming a constant cosmic ray density. Reference was given to 
some of the papers that had used this approach in that section. Somewhat more 
surprising perhaps is the good agreement that is obtained by assuming a 
constant cosmic ray density in the plane for the region (10° < i < 100°, |b| < 
10°) and the COS-B data for the y-ray energy Interval 300 MeV < E < 5000 MeV 
as shown by Lebrun et al. (1983). There are, however, several factors which 
contribute. First, the relatively small percentage of the Compton component 
In the galactic center region compared to higher 1 values partially compen- 
sates for the stronger enhancement that would otherwise occur as 1 decreases 
due to an increased cosmic ray Intensity* Second, in both the work here and 
that of Lebrun et al. (1983) the normalization of the molecular hydrogen 
density is treated as an adjustable parameter. Since the molecular hydrogen 
density Is concentrated towards the galactic center, a large normalization 
value for molecular hydrogen is essentially the same as assuming a positive 
cosmic ray gradient towards the center in terms of the y-rays produced. 

Lebrun et al. (1983) used 3.1 x 10 20 molecules cm" 2 k" 1 $ compared to the 
smaller value used here. Also, assuming an effective area of 47 cm 2 for COS- 
B, their deduced source function for cosmic ray Interactions with matter for 
(300 MeV < E y < 5000 MeV) Is 0*60 x 10" 26 y (H atom)" 1 s" 1 sr" 1 or 7.5 x 10" 26 
q (H atom)” 1 s" 1 , a bit higher than the value of 6.7 x 10" 26 y (H atom)" 1 s" 1 
based on high energy particle physics results, but this difference Is in the 
direction expected since ignoring the Compton radiation has the effect of 
raising the deduced source function. It can also be noted that the variation 
in the observed y-ray Intensity between 50° and 10° In 1 is less for 300 MeV < 
E y < 5000 MeV than in the other two COS-B energy ranges. This feature, If 
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confirmed, may be due to a greater contribution of other sources at lower 
energy toward the direction of the Inner Galaxy* 

The constant cosmic ray assumption combined with the source functions 
based on high energy physics results predicts too Wge a diffuse Y-ray 
Intensity In the anticenter direction by a ratio of about 4:3 to 3:2, as noted 
several times before, e.g., Houston and Wolfendale (1982)* This result Is 
expected, since as the matter density decreases In the outer galaxy* the 
cosmic ray density must also, since there Is not then sufficient gravitational 
attraction to hold the local cosmic ray density. 

i 

This Is an appropriate point to reiterate that there Is also an 
unresolved point source contribution to the "diffuse” radiation measured by 
the SAS-2 and COS-B Y-ray Instruments because of the limited angular 
resolution of these Instruments. It Is quite difficult to estimate this 
contribution; however, several factors suggest that point sources may not be a 
major contributor (see, for example, Cesarsky, 1980). These Include the 
apparent near uniformity of the energy spectrum and the Y-ray luminosity of 
the galaxy and Its distribution being about what would be expected from the 
diffuse sources. For the purpose of this paper, the *~ vader Is simply asked to 
keep In mind that there Is some point source contribution yet to be determined 
which at least for the moment appears to be small. 

IV* CONCLUDING REMARKS 

It was the Intent of this paper to determine If, In light of the recent 
developments, the calculation of the diffuse Y radiation including 
contributions from cosmic ray nucleon-nucleon Interactions, bremsstrahlung, 
and Compton Intensities leads to a reasonable agreement with the diffuse Y-ray 


results when using the current available estimates of the relevant Inter- 
stellar parameters and certain specific assumptions. The two most Important 
assumptions, beyond the acceptance of the Interstellar matter and photon data 
used here and the Interaction cross section Information, are that the relevant 
components of the Interstellar matter all lie In a common spiral pattern and 
that the cosmic ray density Is proportional to the matter density on the scale 
of the spiral arms. These assumptions are supported by observations and 
theoretical considerations as described earlier. In general, the agreement 
between the theoretical predictions and the y-ray data seems reasonable. For 
the present, more detailed refinements are Inappropriate In view of the limita- 
tions Imposed by the data and the limited knowledge related to some of the 
input parameters. 

There are several further conclusions or suggestive Implications which 
emerge. First, since the results obtained do suggest that the general 
concepts are reasonable, as y-ray data of better angular accuracy and energy 
resolution and greater statistical weight become available, It should be 
possible to deduce the galactic cosmic ray and matter density distributions on 
a broad scale and even In relation to clouds with greater accuracy than has 
previously been possible. It may also be possible to detect the arms and even 
specific large clouds on the far side of the galaxy. 

As with the earlier papers noted previously, the variation of the cosmic 
ray density with matter density on a broad arm scale seems consistent with the 
data. The anticenter region seems 'cl early to show a decrease in cosmic ray 
intensity, while the Inner galaxy is consistent with an Increase. At high 
latitudes and mid-longitudes where the cosmic rays are either local or would 
be expected to have characteristics similar to the local cosmic rays, 
agreement between observed y-ray Intensities and those predicted by a uniform 
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cosmic ray distribution seem reasonable. 

The observed energy spectrum Is of some concern. The relatively small 
amount of galactic emission observed by the COS-B Y-ray Instrument above 300 MeV 
places rather significant constraints on models of diffuse galactic emission. 
First, even assuming that essentially all of the galactic y radiation above 
300 MeV not clearly Identified with point sources comes from cosmic ray 
Interactions, a relatively steep cosmic ray electron energy spectrum, only 
marginally consistent -with the spectrun deduced from radio observations, Is 
still required to compensate for the value of the normalization for the 
molecular hydrogen (which Is primarily concentrated In the inner galaxy) 
forced by the high energy Y-ray observations. There are several possible 
explanations for this situation: 

1. Point sources account for the rest of the y radiation In the 
lowest energy Intervals. The few observed Y-ray source spectra are relatively 
hard like the diffuse spectrum, but these hard spectra sources may have been 
observed first for that very reason since the direction of the higher energy 
r-rays can be measured more accurately. If the molecular hydrogen density Is 
even lower, the point sources would become even more Important. If the cosmic 
ray electron spectrum Is closer to the modulated one, as some modulation 
theory would suggest, point sources at lower energies would also be more 
Important. (It might be noted that If this Is the explanation, the Y-ray 
emission spectrum in the outer galaxy* where point sources are presumably much 
less common, would be flatter than observed locally.) 

2. The electron spectrum Is quite steep and intense with the 
molecular hydrogen density being lower. The difficulties with this assumption 
are primarily related to constraints set by cosmic ray and radio data. 

3. The diffuse galactic Y-ray Intensity above 300 MeV derived from 
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the COS-B data Is underestimated for some reason, or the threshold of the hf qh 
energy Interval Is underestimated. It Is only the Intensity above 300 MeV 
which has forced these alternate considerations, and only a small correction 
In the energy threshold would eliminate the difficulty. 

The Compton radiation Is calculated to account for from 7% to 14% of the 
Y radiation In the plane (lb I <10°) and for 8% to over 20% depending on 
energy at higher latitudes. In the central part of the galaxy Its contri- 
bution Is much smaller than It would be If the local electron energy spectral 
shape existed there; In this central region the high energy part of the 
electron spectrum from which the Compton radiation comes has Itself been 
suppressed by the Compton radiation. 

The Compton radiation appears then to account possibly for a major part 
of the approximately 20% difference between the calculation of source 
functions based on high energy physics results, as used here, and the 
observations or source functions calculated by other authors based on the 
simplifying assumption that only cosmic ray matter Interactions are Important 
for the production of the diffuse y radiation. As a result, there is no 
compelling need to assume a large point source contribution (about 10% is 
quite compatible with the work here) although there Is the suggestion that a 
more significant addition may be needed at lower energies than high energies, 
but there might also be an enhanced lower energy cosmic ray electron Intensity 
Hilch would account for this difference. There Is the warning, however, that 
both this paper and the others treated the molecular hydrogen normalization as 
an adjustable parameter and only demanded it be within the fairly wide range 
allowed by other considerations. Even so, the agreement obtained with the y- 
ray data when the molecular hydrogen contribution Is added to the others which 
are not normalized adds credence to the hypothesis that the observed galactic 
diffuse radiation is primarily due to the Interactions of cosmic rays with 


photons and matter. 

We wish to thank the COS-B Caravane collaboration for provldlnq us the 
revised background estimate prior to publication, and Drs. Robert Hartman, 
Bradford Mauger, and Andrew Strong for helpful suggestions. 
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FIGURE CAPTIONS 

Figure 1; The high energy E > 100 MeV Y-ra,y Intensity as a function of 
longitude for -10° < b < 10° from the SAS-2 data (Hartman et 
al., 1979) compared to the model discussed here. 

Figure 2 : Gamma ray Intensity as a function of longitude averaqed over the 

latitude range -10° < b < 10° from 70 MeV - 150 MeV, 150 MeV - 
300 MeV, and 300 MeV - 5000 MeV from the COS-B data (Mayer- 
Hasselwander et al., 1982) compared to the model discussed here 
shown by the solid line. The dashed line In the 70-150 MeV 
graph represents the new "0.0 M line based on the revised 
background Intensity for the COS-B data In this energy Interval 
discussed In the text. 

Figure 3: Energy spectrum of the galactic y radiation for a r&gion near 

the galactic center. The calculated spectra are based on the 
work described here. The solid curves give the sum of all 
components and the two principal components. The dot-dash curve 
Includes an estimated correction for the Increased energy loss 
by electrons In the Inner galaxy. The 300 to 5000 MeV point of 
COS-B (Mayer-Hasselwander et al., 1983), which covers a large 
range In energy Is plotted at an energy where the differential 
energy spectrum of the equivalent power law spectrum Is equal to 
the Integral Intensity divided by the energy Interval width. 

The Compton component shown as a lightly dashed line Is seen to 
be small and uncertain because of the large effect of the 
Compton radiation on the parent electron spectrum in the 
galactic center region as discussed in the text. The COS-B data 
are those of Mayer-Hasselwander et al. (1982), and the SAS-2 
data are those of Hartman et al. (1979). 

Figure 4: Latitude distribution of the diffuse galactic y radiation for 

three energy ranges and three longitude reqions. The data are 
from the experiment (Mayer-Hasselwander et al., 1982); the 
solid lines are the unnormalized predictions of the model 
discussed in the text. The horizontal dashed lines In the 
lowest energy bln are the revised zero level, taking into 
account the more recent estimate of a higher background for the 
energy Interval (Strong, 1982). 
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